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Abstract 

In this paper, we first briefly review some recent results on the distri- 
bution of the maximal eigenvalue of a (A^ x A'^) random matrix drawn from 
Gaussian ensembles. Next we focus on the Gaussian Unitary Ensemble 
(GUE) and by suitably adapting a method of orthogonal polynomials de- 
veloped by Gross and Matytsin in the context of Yang- Mills theory in two 
dimensions, we provide a rather simple derivation of the Tracy- Widom law 
for GUE. Our derivation is based on the elementary asymptotic scaling 
analysis of a pair of coupled nonlinear recursion relations. As an added 
bonus, this method also allows us to compute the precise subleading terms 
describing the right large deviation tail of the maximal eigenvalue distri- 
bution. In the Yang-Mills language, these subleading terms correspond to 
non-perturbative (in 1/N expansion) corrections to the two-dimensional 
partition function in the so called 'weak' coupling regime. 

1 Introduction 

Quite a long time ago, Wigner ]T] introduced random matrices in the context of 
nuclear physics. He suggested that the highly-excited energy levels of complex 
nuclei can locally be well represented by the eigenvalues of a large random 
matrix. A big nucleus is a rather complex system composed of many strongly 
interacting quantum particles and it is practically impossible to describe its 
spectral properties via first principle calculations. The idea of Wigner was to 
model the spectral properties of the complex Hamiltonian of such a big nucleus 
by those of a large random matrix preserving the same symmetry. This was 
a very successful approach in nuclear physics. Since then, the random matrix 
theory (RMT) has gone beyond nuclear physics and has found a wide number 
of applications in various fields of physics and mathematics including quantum 
chaos, disordered systems, string theory and even number theory [2]. A case of 
special interest is the one of Gaussian random matrices (originally introduced by 
Wigner himself) where the entries of the matrix are Gaussian random variables. 
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Depending on the symmetry of the problem, Dyson distinguished three 
classes for the matrix X [3]: 

o the Gaussian Orthogonal Ensemble (GOE) : X is real symmetric, 
o the Gaussian Unitary Ensemble (GUE) : X is complex Hermitian. 
o the Gaussian Symplectic Ensemble (GSE) : X is quaternionic Hermitian. 

Let us write X^ the adjoint of X, i.e. the transpose of X for the GOE, 
the complex conjugate transpose for the GUE and the quaternionic conjugate 
transpose for the GSE. A Gaussian random matrix is a iVx iV self-adjoint matrix 
X, i.e. X'^ — X distributed according to the law 

for GOE 

ViX) cx e-^^(^') with 13^ \ 2 for GUE (1) 

for GSE 

where, for convenience, we have chosen the prefactor /3 of the Tr(X^) to be 
/3 = 1 for the GOE, /3 = 2 for the GUE and /3 = 4 for the GSE. For in- 
stance, for the GUE we have ^ = 2 and V{X) cx e"'^''(^') cx e~^'.^^^'-'^^ as 
X^ = X''X — J2i j l^ij P- This means that X is a, N x N complex Hermitian 
matrix with entries KeXij and ImX^.j for i < j that are independent (real) 
random variables distributed according to the same centered Gaussian law with 
variance 1/4 and the Xi^i are (real) independent Gaussian variables with mean 
and variance 1/2. In case of GSE, there are 2N eigenvalues, each of them 
two-fold degenerate and Tr in ^ for /3 — A is defined so that only one of the 
two fold degenerate eigenvalues in X is counted. 

Self-adjoint matrices can be diagonalized and have real eigenvalues. The 
joint distribution of eigenvalues of the Gaussian ensemble is well known [H [2] 

7'(Ai,...Aw) = Bwe-^^"i^? ni^J-^'^l'^ (2) 

where Bn is a normalization constant such that / (Yii d-K) 'P{Xi, Ajv) = 1 (it 
depends on /3) and the power /3 of the Vandermonde term is called the Dyson 
index /3 = 1,2 or 4 depending on the ensemble (resp. GOE, GUE or GSE). 
Note that we have chosen the prefactor of Tr(X^) term in ([T]) to be the same 
as the Dyson index /3 just for convenience. This prefactor is not very important 
as it can be absorbed by rescaling the matrix entries by a constant factor. In 
contrast, the value of the Dyson index /3 = 1, 2 or 4, characterizing the power 
of the Vandermonde term, plays a crucial role. The normalization constant 
can be computed using Selberg's integral 12^ Bn = — ^ — r(l -|- 

/3/2)^/ [nf=ir(l+j/3/2)]. 

Because of the presence of the Vandermonde determinant Iljxfc i^j ^ ^k) 
in Eq. the eigenvalues are strongly correlated random variables, they repel 
each other. In this paper, our focus is on the statistical properties of the extreme 
(maximal) eigenvalue Amax = max (Ai, A2, . . . , Xn)- Had the Vandermonde term 
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been not there in the joint distribution the joint distribution would factor- 
ize and the eigenvalues would thus be completely independent random vari- 
ables, each with a Gaussian distribution. For such independent and identically 
distributed random variables {A^}, the extreme value statistics is well under- 
stood [5 and the distribution of the maximum, properly shifted and scaled, 
belongs to one of the three universality classes Gumbel, Frechet or Weibull (for 
large N) depending on the tail of the distribution of individual A^'s. However, 
in the case of random matrix theory, the eigenvalues Xi 's are strongly correlated 
variables. For strongly correlated random variables there is no general theory 
for the distribution of the maximum. In case of Gaussian random matrices, 
where the joint distribution ^ is explicitly known, much progress has been 
made in understanding the distribution of Amax following the seminal work by 
Tracy and Widom [6l [7] . This then provides a very useful solvable model for 
the extreme value distribution in a strongly correlated system and hence is of 
special interest. 

Let us first summarize some known properties of the random variable Amax- 
Its average value can be easily obtained from the right edge of the well known 
Wigner semi-circle describing the average density of eigenvalues. For a Gaussian 
random matrix of large size N, the average density of eigenvalues (normalized 
to unity) pn{X) — (-^ i5(A — Ai)) has a semi-circular shape on a finite support 
[—^2N,^2N] called the Wigner semi-circle [1]: 

P7v(A) ~ — = g I — I with g{x) = —Jl — x^ for large iV (3) 

The quantity pN{X)dX represents the average fraction of eigenvalues that lie 
within the small interval [A, A + dX]. Therefore, Eq. ([3]) means that the eigen- 
values of a Gaussian random matrix lie on average within the finite interval 
[—V2N,V2N]. Note also that one can rewrite, using the joint distribution in 

m 

PNiX) = {^Y1 '^(^ - ^^)) = / ^2 . . . , A^) dA2 . . . dXN. (4) 

i 

Hence the average density of states Pn{X) can also be interpreted as the marginal 
distribution of one of the eigenvalues (say the first one). Thus, the marginal 
distribution also has the shape of a semi-circle. Figure [T] shows the average 
density pn{X) (a = 1 here). 

It then follows that the average value of the maximal eigenvalue Amax is 
given for large N by the upper bound of the density support: 

(Amax)«V2iV for large (5) 

However, Amax fluctuates around this average value from one realization to 
another and has a distribution around its mean value ^/2N (see Fig. [1] with 
a — I). What is the full probability distribution of Amax? From the joint 
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Figure 1: Average density of the eigenvalues of a Gaussian random matrix 
Pat (A) as a function of A (blue dashed line). The density has a semi-circular 

shape ( "Wigner semi-circle" ) and a finite support [— ^J^^] ■ The maximal 

eigenvalue has mean value (Amax) ~ \J f'^'" large N and its distribution close 

to the mean value, over a scale of 0(A^^^/^) has the Tracy- Widom form (red 
solid line). However, over a scale (V^iV) the distribution has large deviation tails 
shown by solid green (left large deviations) and solid blue (right large deviations) 
lines. 



distribution of eigenvalues in Eq. ([2]), it is easy to write down formally the 
cumulative distribution function (cdf) of Amax as a multiple integral 

(Amax < i) = n / '^^^ n i^j- - ^fei'' (6) 

which can be interpreted as a partition function of a Coulomb gas in presence 
of a hard wall at the location t (see the discussion in Section 2). The question is 
how does Vn (Amax < t) behave for large A^? It turns out that the fluctuations 
of Amax around its mean yj2N have two scales for large N . While typical 
fluctuations scale as iV~^/^, large fluctuations scale as N^^"^ and their probability 
distributions are described by different functional forms (see Fig. [T]with a — 1). 

Typical fluctuations: From an asymptotic analysis of the mutilple integral in 
Eq. ([6]) , Forrester j8j , followed Tracy and Widom [6l [7] deduced that for large 
N ^ small and typical fluctuations of the maximal eigenvalue around its mean 
value \/2N are of order 0(iV~^/^) and can be written as 

Amax~\/2A-f a;37V-Sx (7) 

where ai,2 = l/%/2 (for GOE and GUE) and 04 = 2-'^/^ (GSE) and x is a 
random variable characterizing the typical fluctuations. Tracy and Widom |6l[7] 
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proved that for large A'', the distribution of x is independent of A'': P (x ^ a;) = 
Fi^{x). The function i^^(a:) depends exphcitly on /? and is called the Tracy- 
Widom distribution. For example, for /3 = 2 [7] , 



^2(2;) = exp 



[z — x)q^{z)dz 



(8) 



where q{z) satisfies the special case of a = of the Painleve II equation 

q"{z) = 2q^{z)^ zq{z) + a. (9) 

For a = 0, the solution only requires the right tail boundary condition for its 
unique specification: q{z) ^ Ai(z) as z — >■ 00, where Ai(z) is the Airy function 
that satisfies the differential equation Ai"(z) — zAi(z) — and vanishes as, 
Ai(z) « 2^z^/'^ e~^^^'^ as z — >• 00. This solution of the special case a = of 
the Painleve-II equation is called the Hastings-McLeod solution [S]. For /3 = 2 
and ;5 = 4, one has [51 [7] 

Fi{x) = [F2{x)]''^ cxp 
F^{x) = [F2{x)]^''^ cosh 



q{z)dz 
q{z)dz 



(10) 
(11) 



Note that Fp{x) — Prob(x < x) is the cumulative probability of the scaled 
random variable x ^-nd hence it approaches to 1 as a; — ^ cx) and vanishes to 
as X — —00. The corresponding probability density function (pdf) F'^^x) = 
dFjs (x) I dx vanishes as a; — >■ ±00 in an asymmetric fashion 



F'M 



exp 



exp 



24 ' ' 



-CX) 



2/3 



,3/2 



as x 



(12) 
(13) 



Over the last decade or so, the Tracy- Widom distribution has appeared in a 
wide variety of problems ranging from statistical physics and probability theory 
all the way to growth models and biological sequence matching problems (for 
reviews see [lOl [TTl [121 US])- These include the longest increasing subsequence 
or the Ulam problem [lH [151 IlH], a wide variety of (l-l-l)-dimensional growth 
models [16l [17l [HI [TH [20] , directed polymer in random medium |21] and the 
continuum Kardar-Parisi-Zhang equation [22l [231 EH [21], Bernoulli matching 
problem between two random sequences [2 6) , nonintersecting Brownian motions 
(see e.g. [27l[28] and references therein). This distribution has also been mea- 
sured in a variety of recent experiments, e.g., in the height distribution of fronts 
generated in paper burning experiment [29j , in turbulent liquid crystals 30^ and 
more recently in coupled fiber laser systems j31) . 

Large deviations: Tracy- Widom distribution describes the probability of typ- 
ical fluctuations of Amax around its mean (on a scale of but not the 
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atypical large fluctuations, i.e., fluctuations of order 0{vN) around the mean 
value \/2N. Questions regarding such large/rare fluctuations do arise in vari- 
ous contexts [32j |33l |34] and have recently been computed [32l [33l [Ml [35] to 
dominant order for large N. As a summary, the probability density of Amax, 
V (Amax = = fti^N (Amax < t)], IS given for large N by: 



t) 



exp 



for t < V2N and \t 



/2N\ w 0{VN) 



pf I t-V2N 



($^) for 1^ - V2]V| « O (^-Va) 



exp |-/3A^V'+ (^) + •■•} for t > V2N and \t 



/2N\ 



(14) 

where Fp (x) is the Tracy- Widom distribution and where ip- and "0+ are respec- 
tively the left and right large deviation functions describing the tails of the dis- 
tribution of Amax- The rate function i}-{z) was explicitly computed in [32ll33] . 
while ip+i^) was computed in |35) . both by simple physical methods exploiting 
the Coulomb gas analogy. A more complicated, albeit mathematically rigorous, 
derivation of ip+{z) in the context of spin glass models can be found in [55] . 
These rate functions read 



V'+(^) 



z 

y 



108 



[z-^ + 15z) 



108 



z\fz^ 



+ ln 



V2 



+ 6 + z 



V2 



for z > V2. 



In 3 



, for z < a/2 



(15) 



Note that in ref. 135' , the function '0+(z) was expressed in terms of a com- 
plicated hypergeometric function, which however can be reduced to a simple 
algebraic function as presented above in Eq. p5]) . Note also that while Fp{x) 
depends explicitly on /3, the rate functions ij^-iz) and V'+(^) are independent 
of (3. These rate functions only give the dominant order for large N in the 
exponential. In other words, the precise meaning of « is that for large N: 

1™ UTTtIhT' (Amax = z^/n] = -V-(2)forz < \/2and lim -^-rrlnT' (Amax = z^/N] = ~il)+{z) 

for z > \f2. When z approaches -\/2 (from below or above) it is easy to see that 
the rate functions vanish respectively as 



i>-{z) -> 



1 



6^2 

27/4 



2-z)^ + ... 



V^+(z) ^ _(z-V2)3/2 + 



(16) 
(17) 



Note that the physics of the left tail [32] [33] is very different from the physics 
of the right tail [35]. In the former case, the semi-circular charge density of the 
Coulomb gas is pushed by the hard wall {z < \/2) leading to a reorganization 



6 



of all the N charges that gives rise to an energy difference of 0{N^) [351 155] . 
In contrast, for the right tail z > v^, the dominant fluctuations are caused by 
pulling a single charge away (to the right) from the Wigner sea leading to an 
energy difference of 0{N) [55] , 

The different behaviour of the probability distribution for z < \/2 and z > 
y/2 leads to a 'phase transition' strictly in the iV — )■ oo limit at the critical point 
z = \f2 in the following sense. Indeed, if one scales A by \/]V and takes the 
N ^ oo limit, one obtains 

- lim lnPjv(A„,ax < ^ViV) = ^^(z) for z < V2 

= for z > %/2 (18) 

Note that since Fjv(Aniax < t) can be interpreted as a partition function of a 
Coulomb gas (see Eq. ([6])), its logarithm has the interpretation of a free energy. 
Since i^-{z) ~ {\/2 — z^ as z — )■ ^/2 from below, the 3-rd derivative of the 
free energy is discontinuous at the critical point z = \/2. Hence, this can be 
interpreted as a third order phase transition. 

However, for finite but large N, it follows from (|14|) that the behavior to the 
left of z = \/2 smoothly crosses over to the behaviour to the right as one varies 
z through its critical point z = \/2 and the Tracy- Widom distribution in (|14l) 
around the critical point is precisely this crossover function. Indeed, if one zooms 
in close to the mean value y/2N by setting t — \/2N +xN^'s / (for /3 = 1, 2) in 
the rate functions ^^{t/\/lV) and ip+lt/^/N) in ([T4|) . one expects to recover, by 
taking large N limit, respectively the left and the right tail of the Tracy- Widom 



3/2 



3N 



distribution. With this scaling, and using pT)) . one finds V'+(i/v^) 
and thus V (A^ax = i) ~ exp | — ^x"^/^!, which indeed matches the dominant 
order in the far right tail of the Tracy- Widom distribution for /3 = 1 , 2 in (jTF 



Similarly for the left tail {x < 0), using (|16p . one finds ip^(t/\/N) 2lW^ thns 

V (Amax = t) ^ exp | — ^l^^l'^l which matches the left tail of the Tracy- Widom 
distribution in ((T^ . 

More recently, higher order corrections for large N have been computed for 
the left tail of the distribution jiJI using methods developed in the context of 
matrix models. Note that in [37] a different notation for /3 was used: /3 = 1/2 
(GOE), 13 = 1 (GUE) and /3 = 2 (GSE). To avoid confusion, we present below 
the results in terms of the standard Dyson index /3 = 1, 2, 4. 

V (Amax = i) « exp |-$A, (^-^,13^ I for t < V2N and \t - V2N\ w 0{Vn) 

(19) 

where 

$w(z,/3) = /3iV2^_(z) + 7V(/3-2)$i(z) + 0^ lniV + $o(/3,^) (20) 
with ip-iz) given in Eq. ()15p (dominant order). The subleading terms are given 
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by m 



$i(z) = + _(z2 + 6)4(7^2 + 6-2z)2 

In 18 1 , 



and 



2^/^^^ - z + V3{z^ + 6)i - 2z)^ (21) 

(22) 



and (see Eq. (4-35) in Ref. ,37 
1 „ 1 1 



*o(/3,^) = 



- , , o / 19 19/3 9\ , „ 

/3H n2+ \ + - n3 

12^ 3^ 'Sj \12I3 48 8/ 
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3 J_ ^ 

8 ^ 4^ ^ 16 
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2 3^ 12 



V6^ 



3^ 3 



In 



\l-2z + v/eTI^ + (z^ + 6) 



3 V /3 4 ' 



-2z + Ve + z^ 



-z + 2^6 + z2 + V3 + 6) ^^■^ \/ V6 + z2 _ 2^ 



- In 



(-18 + z2) z+ (6 + z2) 



2n3/2 



ln/3 



where np is a complicated function of (3. For (3/2 integer, it reduces to [37] 



(23) 



/^/2 + l\ 2C^ ln(/3/2) ^ 2m / m 



For instance, for the GUE (/3 = 2), we find K2 = + C'(-l)- For /3 = 1, 2 

and 4, the expression in Eq. (|20|) matches the left asymptotics of the Tracy- 
Widom distribution, i.e. the asymptotic behaviour of Fl^{x) for a; — >■ — oo, see 
ref. [35] • However, for the right tail of the distribution of Amax, the corrections 
to dominant order for large N are, to our knowledge, not known until now. In 
fact, one of the results of this paper is to compute these right tail corrections 
for the GUE {/3 — 2). Both left and right large deviations are plotted in Fig. [2] 
for the GUE. The left tail is described by $Ar(z, 2) in Eq. jSO]), the right tail is 
described by our result given in Eq. (126^ . 

Another result of this paper concerns a simpler and pedestrian derivation of 
the Tracy- Widom distribution for the GUE case. The original derivation of the 
Tracy- Widom law for the distribution of typical fluctuations of Amax [6] [7] is 
somewhat complex as it requires a rather sophisticated and nontrivial asymp- 
totic analysis of the Fredholm determinant involving Airy Kernel [S] [7] . Since 



(24) 
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this distribution appears in so many different contexts, it is quite natural to 
ask if there is any other simpler (more elementary) derivation of the Tracy- 
Widom distribution. In this paper, we provide such a derivation for the GUE 
case. Our method is based on a suitable modification of a technique of or- 
thogonal polynomials developed by Gross and Matytsin [33] in the context of 
two-dimensional Yang-Mills theory. In fact, the partition function of the con- 
tinuum two-dimensional pure Yang-Mills theory on a sphere (with gauge group 
U(A'^)) can be written (up to a prefactor) as a discrete multiple sum over inte- 
gers [101 m] 

oo 

Z{A,N)= n {n,-n,)'e-^^/^^^^"-^''" (25) 

Tii,n2,...,njv = — oo l<i<j<N 

where A is the area of the sphere. In the N oo limit, the free energy InZ, as 
a function of A, undergoes a 3rd order phase transition known as the Douglas- 
Kazakov transition [32] at the critical value Ac — tt^. For A > Ac, the system 
is in the 'strong' coupling phase while for A < Ac, it is in the 'weak' coupling 
phase. For finite but large N, there is a crossover between the two phases as one 
passes through the vicinity of the critical point. In the so called double scaling 
limit (where A — >• Ac, iV — >• oo but keeping the product (A — Ac)N'^^^ fixed), 
the singular part of the free energy satisfies a Painleve II equation [35] • Gross 
and Matytsin (see also [13]) used a method based on orthogonal polynomials 
to analyse the partition sum in the double scaling limit, as well as in the weak 
coupling regime {A < Ac) where they were able to compute non-perturbative 
(in 1 /N expansion) corrections to the free energy. Actually, a similar 3-rd order 
phase transition from a weak to strong coupling phase in the N ^ oo limit 
was originally noticed in the lattice formulation (with Wilson action) of the two 
dimensional U(iV) gauge theory [44j|45l|T4] and in the vicinity of the transition 
point the singular part of the free energy was shown to satisfy a Painleve II 
equation [15]. Note that similar calculations involving the asymptotic analysis 
of partition functions using orthogonal polynomials were used extensively in the 
early 90 's to study the double scaling limit of the so called one- matrix model 
(for a recent review and developments, see e.g. Ref. [47J. 

In our case, for the distribution of Amax, we need to analyse the asymptotic 
large N behaviour of a multiple indefinite integral in Eq. ([6]) , as opposed to the 
discrete sum in Eq. ([25]). However, we show that one can suitably modify the 
orthogonal polynomial method of Gross and Matytsin to analyse the multiple 
integral in Eq. ([5]) in the limit of large N. In fact, we find a similar third order 
phase transition (in the N ^ oo limit) in the largest eigenvalue distribution 

(Amax < i) as a function of t at the critical point tc ~ V27V. The regime of 
left large deviation of Pat (A,nax < i) {t < tc) is similar to the 'strong' coupling 
regime {A > Ac) of the Yang-Mills theory, while the right large deviation tail 
of Pjv (Aniax < t) {t > tc) IS similar to the 'weak' coupling regime {A < Ac) of 
the Yang-Mills theory. For finite but large N, the crossover function across the 
critical point that connects the left and right large deviation tails is precisely 
the Tracy- Widom distribution. Thus the Tracy- Widom distribution corresponds 
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precisely to the double scaling limit of the Yang-Mills theory and one finds the 
same Painleve II equation. A similar 3rd order phase transition was also found 
recently in a model of non-intersecting Brownian motions by establishing an 
exact correspondence between the reunion probability in the Brownian motion 
model and the partition function in the 2-d Yang-Mills theory on a sphere [371 

SB]. 

The advantage of this orthogonal polynomial method to analyse the max- 
imum eigenvalue distribution is twofold: (i) one gets the Tracy- Widom distri- 
bution in a simple elementary way (basically one carries out a scaling analysis 
of a pair of nonlinear recursion relations near the critical point and shows that 
the scaling function satisfies a Painleve II differential equation) and (ii) as an 
added bonus, we also obtain precise subleading corrections to the leading right 
large deviation tail (t > V^N). The subleading corrections, in the Yang-Mills 
language, correspond to the non-perturbative corrections in the weak coupling 
regime as derived by Gross and Matytsin [39]. More precisely we show that 



r (Aniax - i) ~ e~^^^+^Vw) for t > V2N , t - V2N « 0{V 

(26) 

where ^+{z) is given in Eq. (|15p . Note that only the leading behaviour 
exp[—2N'ip+{z)] was known before [35], but the subleading corrections are, to 
our knowledge, new results. We also verify that our expression matches the 
precise right asymptotics of the Tracy- Widom distribution. Figure [2] shows the 
distribution of Amax for the GUE: close to the mean value it is described by the 
Tracy- Widom distribution, whereas the tails are described by the large devia- 
tions. The right tail (right large deviation) is given by our result in Eq. (|26|) . 
Together with the subleading terms in the left tail in Eq. (|20|) . our new result 
in Eq. I^E^ then provides a rather complete picture of the tail behaviors of the 
distribution of A,nax on both sides of the mean ^/2N . 

The rest of the paper is organized as follows. In Section [21 we start with 
some general notations and scaling remarks for the GUE. In Section [3l we 
explain the method of orthogonal polynomials on a semi-infinite interval and 
derive some basic recursion relations. In Section [31 we compute the right tail 
of the distribution of Amax (dominant order and corrections for the GUE): it 
describes atypical large fluctuations of Amax to the right of its mean value. In 
Section [5l using results of the previous sections and basic scaling remarks, we 
derive the Tracy- Widom law (with /? = 2 for the GUE) that describes small 
typical fluctuations close to the mean value. 



2 Notations and scaling 

In the rest of the paper we focus only on Gaussian random matrices X drawn 
from the GUE (/3 = 2). They are Hermitian random matrices — X such that 

ViX) (X e""^^^^ ) where we have introduced an additional parameter a > 
for the purpose of certain mathematical manipulations that will be clear later. 
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Figure 2: Rate function — ln7'(Amax = za//V) associated to the distribution 
'P(Amax = of the maximal eigenvalue of a random matrix from the GUE for 
large N . Close to the mean value z = v^, the distribution is a Tracy- Widom 
law (red line), it describes the small typical fluctuations around the mean value. 
Atypical large flutuations are described by the large deviations: the left large 
deviation in green (z < v^), the right deviation in blue {z > V^). 
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Setting Q = 1 at the end of the calculation, we will recover the usual GUE. 
With the additional parameter a, the joint distribution of the eigenvalues of X 
is given by (see Eq. 

V{Xi,...,XN)^BN{a)e-"^^^^" UiXj-Xkf (27) 

j<k 



is the normalization constant. The 



where BAr(a) = (2a) (27r)-f / Il^^i n\ 
Vandermonde determinant appears with a power 2 as we consider the GUE 
{(3 ~ 2). This determinant indeed comes from a Jacobian for the change of 
variables from the entries of the matrix to its eigenvalues. The power is related 
to the number of real degrees of freedom of an element of the matrix, which is 2 
for complex entries, i.e., for GUE (it is 1 for real entries GOE and 4 for quater- 
nion entries GSE). As we will see later, this power 2 is crucial for the method 
of orthogonal polynomials to work. The technique of Gross and Matytsin 
that we adapt here also works only for the GUE /3 = 2 case. 

Given the joint distribution of eigenvalues in Eq. ([27]), it is easy to write 
down the cumulative distribution of the maximal eigenvalue Amax 

Pat (A„,ax < y) = Prob [Xi < y, X2 < y, . . . , Xn < y] = ^lS^I^ (28) 

where the partition function Z^r is given by the multiple indefinite integral 

The normalization Zjv(oo,ck) is actually related to BN{a) in Eq. (1271) as 
Bisf{a) = 1/ {Nl Zpf{oo,a)). Note that, by the trivial rescaling y^Xi — > 
in dMl), it follows from (EHl) that 



TD> n ^ \ ^N{y,a) ZN{yy/a,l) ,on^ 
Pat (Amax < y) = :^-7 r = -^-^ r^- (30) 

Thus y and a always appear in a single scaling combination y\foL. 

We will henceforth focus on the large iV limit. For fixed a, one can easily 
figure out from the joint pdf in Eq. ([27| how a typical eigenvalue Atyp scales 
with 'N for large N . Let us rewrite the joint distribution of eigenvalues in Eq. 
dUl) as 

P(Ai,...A^)«expJ -a^A2 + 2^1n|Aj-Afe| I (31) 

which can then be interpreted as a Boltzmann weight oc exp[— i?eff]) with ef- 
fective energy E^s = a ^ ■ Af — 2 J2j<k l-^i ^ ^k\- The eigenvalues can thus 
be seen as the positions of N charges of a 2D Coulomb gas (but restricted to 
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be on the real line) which repel each other via a logarithmic Coulomb potential 
(coming from the Vandermonde determinant in Eq. ((27)) ) [T2] . In addition, the 
charges are subjected to an external confining parabolic potential. For large N, 
the first term of the energy is of order iVA^yp, whereas the second is of order 
A^^ because of the double sum. Balancing the two terms NX^yp ~ 7V^ gives the 
scaling of a typical eigenvalue for large N: Atyp ~ VN- For large iV, the eigen- 
values are close to each other and they can be described by a continuous charge 
density (normalized to unity) pn{X) = j; S {X — Xi) . The average density 
of states for large N can be obtained by evaluating the full partition function 
Zn{oo, a) (the denominator in Eq. (^5]) ) via a saddle point method. The saddle 
point density is the density that minimizes the effective energy (see the book 
of Mehta [1]) E^s = aN J dXpN{X)X^ - J dX J dX' pNiX)pNiy)\n\X - X'\ 
(in its continuous version). This gives the well-known semi-circle law (which is 
exactly the same as in Eq. ^ for a = 1): 



PnW-^p(-^) withp(.) = ^y^-.^ (32) 



The density is plotted in figure [TJ 

The average value of Amax is again given for large N by the upper bound of 
the density support (see Fig. [1]): 

(A,„ax)«\/ for large TV (33) 

V a 

For a = 1, this evidently reduces to the usual expression for (Amax)- The typical 
scaling for large N is thus Amax ~ V^- Hence, we will use Amax = z Vn, where 
z is of order one. 



3 Orthogonal polynomials 

In this section, we introduce the method of orthogonal polynomials to eval- 
uate the partition function in Eq. (1291) . As mentioned in the introduction, 
to evaluate this multiple indefinite integral we will adapt the method devel- 
oped by Gross and Matytsin [35] to enumerate the partition sum (P5|) in the 
two-dimensional Yang-Mills theory. Orthogonal polynomials are very useful to 
handle the square Vandermonde determinant in the distribution of the eigen- 
values in Eq. (j27p . A Vandermonde determinant can indeed be written as 

Ili<ji^j ~ ^i) = det (^r^) . = det (pj-i(Ai)) where pj{X) = X^ + ... is any 

polynomial of degree j with leading coefficient one. The idea is to choose well 
these polynomials pj in order to simplify the computation of the integral. 
We define an operation on pairs of polynomials as follows: 

(/,g>= r dAe-"^' /(A)5(A) (34) 

'J —oo 
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We consider a family {pn}„>o orthogonal polynomials with respect to the 
operation defined above, i.e. with weight e"""^ on the interval ]—cx), ?/]. Without 
any loss of generalization, we define the polynomial PnW of degree n such that 
the coefficient of A" term is always fixed to be 1, i.e., J3„(A) = A" + .... These 
polynomials satisfy the orthogonality property: (pmPm) — for all n ^ m. We 
also write /i„ — {pn,Pn)- Thus, 

{Pn,Pm) = Sn,m K for all n > (35) 

Note that Pn(A) and /i„ are implicitly functions of a and y, i.e. Pn(A) = 
Pn{X\y,a) and /i„ = hn{y,a). In particular, we can easily compute by hand 
the first few p„(A)'s for fixed a > and y, but the expressions become rather 
complex as n increases and it is hard to find a closed form expression for p„ (A) 
for every n (except for the limiting case ?/ — >■ oo): 



Po(A) = 1 
Pi{\) ^ A + 



-ay 



/na [1 + erf (y-^a)] 



-2yV^-e^°0F(l + 2y^a)(l+crf[yy^) 

P2(A) = A H + 

^/^{2 + 2e!''"0FyV" (1 + erf [yV"]) - e'^y^°'TT (1 + erf [yVo]) 



^ -4 - 4e^'"0Fy yg (1 + erf [y^) + e^^'°7r (1 + erf [y^/^)^ 
4a + 4e2^''"V^ya3/2 (1 + erf [y^]) - 2e2y^"7ra (1 + erf [y^])^ 

Thus we get for instance ho ~ (poIpo) = [1 + erf(yv^)] 

and hi = {pi\pi) = -^^2 v^fi+"f(«v^)l ^ etc. In the limit 

y — >■ c», we recover the Hermite polynomials: po ~ I, pi ^ X, P2 = ~ and 

3.1 Partition function 

The partition function ZN{y,a) in Eq. ()29p can be written as a function of 
the hnS. By combination of rows, the Vandermonde determinant in the joint 
distribution of the eigenvalues can indeed be written 



n(Afe - A,) = det (Ar'j^ . = det (p,-i(A,)) 
j<k 
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Then, the partition function can be expressed as 



1 ^ r 

— Jly dA, det(pj_i(A,))^^^. det(p;_i(Afe));^._; e" 
dXe'"^" K„i(A)p,_i(A) 

= det = det ((5ij/ii_i)j = 



det 



i=0 



where in going from the second to the third Une we have used the Cauchy-Binet 
formula [2]- Note that this step works only for /3 = 2. Therefore the partition 
function reduces to: 



ZN{y,a) = Y]_ ^n{y,a) 



(36) 



Thus the integral ZN{y, ce) is now expressed as a product of the coefficients 
The goal of next subsection is to find recursion relations for the /i„'s in order 
to compute them and subsequently analyse their product Zj^ in Eq. in the 
large N limit. 



3.2 Recursion relations 

In general, for orthogonal polynomials (with any reasonable weight function), 
one can write a recursion relation of the form; 

Ap„(A) = p„+i(A) + SnPnW + i?„p„-i(A) (37) 

where Sn and i?„ are real coefficients. This relation comes from the fact that 
Pn = A" + ... and that (pnlo) = for any polynomial q{X) of degree strictly 
smaller than n. The coefficients 5„ and i?„ are functions of a and y, i.e. 
Sn = Sniy, a) and i?„ = R^iy, a). 

Let us first demonstrate that the coefficients i?„ and Sn can be expressed in 
terms of /i„'s. From Eq. ([371), ^e get: (p„_i|Ap„) = i?„ {pn-i\Pn-i) = Rn hn-i- 
On the other hand, we have (p„_i|Ap„) = (Ap„_i|p„) = (p„ + Sn-iPn-i + 
Rn-iPn-2\pn) = {pn\pn) = hn- Therefore RnK-i = hn, thus i?„ = hn/hn-l- 

From Eq. ([31]) again, we also get: (p„|Ap„) = Sn {pn\pn) = 5'n/iri. By 
definition we have (p„|Ap„) = /^_^ dAe""^ Ap^j (A). After integrating by part 
we find: (p„|Ap„) = -^Lg-^J^^pS ^ ..L^^h^. The last equality fol- 
lows from the definition of As /i„ = j^^d\e~°'^ p\ (A), we have indeed 
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= (y|y,a) + 2(p„|%^). However, (p„|^) = as % is a 

polynomial of degree strictly smaller than n (as Pn {X\y,a) = A" + ...). There- 

J_ d In h, 

2a dy 



fore Sn--— 



Combining these relations, i?„ and 5„ are then given by: 

Rn[y,a) = T r and Sn{y,a) = -7; ^ (38) 

hn-i(y,a) 2a dy 

Iterating the recursion relation hn — 1 starting from n — 1, we can write 

hn in terms of i?fc's: 

hn = (^f[ Rk^ ho (39) 

Substituting this result in Eq. we see that the partition function Zn can 

be entirely expressed in terms of a product over the i?„'s. Thus, if we can de- 
termine i?„'s, we can evaluate the partition function explicitly. 



Thus the next task is to determine the i?„'s. To do this, we will first de- 

2fhL = 

da 



rive a set of coupled recursion relations for R^s and ^^'s. We have 



= - /^^ dX e-"^' X^pI (A) ^ -(Ap„|Ap„) where we have used the fact 

that (Pnl^f) = (which follows from the observation that is a polyno- 
mial of degree strictly less than n and hence orthogonality dictates that it is 
identically zero). On the other hand, using Eq. ([57]) we find: (Aj5„|Ap„) = 

(j)n+l + SnPn + RnPn-l\Pn+l + S^Pn + RnPn-l) = h^+l + S'^hn + R^hn-l = 

hn {Rn+1 + S^ + Rn). Therefore, 

- - Rn + Rn+1 + Si 40 

oa 

We can eliminate hn from the relations (p8| and pO)) and get a pair of coupled 
nonlinear recursion relations for i?„ and 5„. Using Eq. (UHl) for n and n — 1, 
and as i?„ = /i„//i„_i, we find - ^'g^" = i?„+i - i?„_i + Sf-^ - Sl_^. Using Eq. 
dSHl), we also find ^^f^ = 2a(5'„_i - S*,!). Finally, we then get our desired 
recursion relations: 

Rn+1 — ^ Rn-i ~ + S^-i (41) 

Sn = Sn-1-^'-^ (42) 

2a ay 

It is easy to show by induction that the two relations (HT|) and (|42p together 
with the initial conditions i?Oj ^1 a-nd Sq uniquely determine i?„ and 5„. The 
additional initial condition Hq is enough to determine /i„ as given in Eq. p9p . 

By definition, po is a polynomial of degree with dominant coefficient 
1. Thus po(A|y,a) = 1. Therefore hoiy,a) = (po|po) = /^^ dA g-"^' = 
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^(1 + erf(?/i/a)). We also have Ro{y,a) = as the recursion relation in 



1 

2 

Eq. ([37| must reduce for n = to Xpo{X) = pi{X) + Sopo{X), i.e., pi{X) 
X - So. Moreover we get from Eq. (|3H1) 5*0 = = ^e-^'y^ /{2aho) = 

—Q-°'V j [y^Tra (1 + erf(yY^))]. We now have i?o, <S'o and ft,o, we can thus deter- 
mine i?i from Eq. (001) for n = 0: i?i = - ~Rq-SI = ySa + l/(2a) - 5*2. 
Thus, the initial conditions can be summarized as: 





= 0, 


Po{X\y,a) 


= 1 , 


ho{y,a) 


J — OO 




1 HF 
2\Ja 


So{y,a) 


1 


9 In ho 






dy 


2aho 


Ri{y,a) 


= ySo - 


^2a ^0 





-ay 



/Tra (1 + ei:i{y^/a)) 

(43) 



3.3 Normalization: limit y — )■ oo 



As y — oo, we can explicitly compute the functions Sm Rn and As men- 
tioned above, in the limit y — oo the polynomials Pn are indeed simply the 
Hermite polynomials. Hence everything can then be computed explicitly in this 
case. We have /io(oo,a) = J^^dX e~"^ = and 5*0(00,0) — 0. Then, by 
recurrence it is easy to show that: 

Sni<x,a) = Rni^.a) = — (44) 

za 

Finally, using Eq. ([39]) we get 



and thus (see Eq. (jSH 



AT 



ZAr(oo, a) = (27r)^ (2a)"^ J| T{n) (46) 

which could also have been computed directly using Selberg's integral. We 
recover the normalization Bj^ = 1/ {Nl Zm) in Eq. (P7|) . 

4 Right tail of the distribution of Amax* large 
deviation function 



In the previous section we derived a pair of coupled recursion equations (j4ip 
and ((42l) with initial conditions given in Eq. p3)) that determines uniquely 
and thus ft,„ and subsequently Zjv via p6p . However, these equations are 
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hard to solve explicitly for general n and y (apart from the case y = +00 as 
shown in the previous subsection). In this section, we derive an approximate 
solution for and hence for the cdf ^E\\ . in the large N limit where N is the 
number of eigenvalues of the matrix X and we will see that this solution for cdf 

is valid only for Amax > (Amax) = \/^' i-^-' "^^^^ describes the right tail of 
the probability distribution. 

We have seen in Section 2 that for large N and fixed a, the maximal eigen- 
value typically scales as Amax ~ V^- We are going to work on this scale, hence 
in the definition of the maximal eigenvalue cdf in ([28l) and p9|) . we will set 
y = zy/N with z w 0(1). We will then work in the limit of large A'' with z fixed. 
With this scaling, the operation (/, g) defined in Eq. for polynomials de- 
pends on N (since the upper limit of integration in (p4| is now z^/N). The 
coefficients i?„, Sn and for a given n, are also now implicitly functions of N. 
We can make an expansion of these parameters for large N and fixed n. The 
dominant order will be given by the y = +00 case (as in previous subsection) 
as y = z^/N 00 as iV ^ 00 for fixed z. In this section, we want to determine 
the first correction to this dominant term. 

However the partition function Zn and thus the cumulative distribution 
(cdf) of Amax is a product of all the hnS ior < n < N. Our expansion will 
thus give us the behaviour of the cdf of Amax for large N only if we can show that 
it is valid not only for fixed n but also for n of order up to N. This constraint of 
validity will be discussed later. We will see that this expansion is actually valid 

only on the right of the mean value, i.e. for y > \J^^ or equivalently z > \J^- 
This method allows us to describe the right tail of the large deviation of the 
distribution of Amax, i-e. 'P(Amax = t) with t > and t - « 0{VN). 

4.1 Expansion of i?„ and Sn 

Let us start by expanding the initial conditions for large N. With the scaling 
y — z\/N with z w 0(1), equations (|43|) become (for z > 0) for large N: 

ho{zVN,a) « .[^ L^e-^"^V... 

1 ~Naz^ 



R,izVN,a) « -l_-£^e-^"^^ (47) 
2a iJaix 



The dominant term for large N corresponds to the limit y — > cx) (see previous 
section): = dXe^"^ = /io(oo,a). Therefore, ignoring the exponen- 
tially small correction for large N leads to Rn{zVN,a) ~ Rn{oo,a) = ^ and 
Sn{zVN,a)^0. 

We want to determine the first correction for large N. Let us make the 
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following ansatz: 



RnizVN, a) « — + c„ e"^""^ and Sn{zVN, a) « d„ e" 
2a 



■Nc 



(48) 



where c„ = Cn{zvN,a) and d„ = dn{zvN,a) are expected to be polynomials 
of z^/N. This will be valid as long as Cn{y,a)e~^°'^ <^ ^ where y = z\fN . 
The initial conditions in Eq. (|Tfl) give: 



co(y,a) = , ci(y,Q:) 



2/ 



2x/Q!7r 



and do{y,a) 



1 



(49) 



Let us replace i?„ and Sn in the recursion equation (pij) and (H^ by the 
ansatz in Eq. (|48l) . We see that and 5,^^^]^ are actually negligible in the 
equation (HT|) giving -R„+i, as they are exponentially smaller than the i?fe's. We 
thus get recursion relations for the coefficients Cn{y,a) and dn{y,a)'. 



2 r 9(ac„) 2 

- i K ay c„ 

n oa 

2ayc„ - — - 
oy J n 



(50) 



h„ is the product of the Rk, it can thus be expressed in terms of Ck- As Zn is 
a function of the /i„'s, it is thus a function of ho and the i?„'s for 1 < n < A^, 
we can use the ansatz (|^ only if c„(za/]V, a)e 
this case we can write: 



Naz ^ n fQj. j^ll n < N. In 

2ot 



lTihn{y,a) = In/io + ^ lni?fc w In/io + ^ 

1 



k=l 



k=l 



2acfc 



-Nc 



ln/i„(oo, a) + 



where y = zV N, and In Z]^{y, a) = Y^^=q hn{y, a) is thus given by 

7V-1 ji 



E 

k=l 



2a 

2ack{y,a) 



(51) 



In Zisi{y, a) « In Zisi(oo, a) + 



TV 



22/ 



EE 

ji=0 fc=l 



2acfc(j/,a) 



-Nc 



(52) 



The partition function only depends on the c^'s. We want now to solve the 
recursion relation for the c^'s in Eq. ([50]) . We do not need to determine the 

dfe's. 



4.2 Solution of the recursion for the c„ 

Let us define ^ and Gn such that: 

^ = ay^ — Naz^ and c„(y,a) 



2V< 



G„(0 



(53) 
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^ is large for large N (proportional to N). G„(^) depends only on f = ay"^. 
This can easily be shown by recurrence with initial condition Gi(^) ~ 1 (as ci 
is given by Eq. (05])). The recursion ([50)1 for the c„'s becomes 

G„+i(0 - G„_i(0 = - I - G„(e) + ^ GUol (54) 



with initial condition Go(C) =0 and Gi(^) = 1. By recurrence again, it is easy 
to show that G„(^) is a polynomial of ^ of degree (n — 1), with leading coefficient 

Let us consider the generating function of the {G„(^)}: 

oo 

F{^,x) = J2^"Gn{0 (55) 

n=l 

The G„(^) are obtained from by a contour integration: 

where G is a contour in the complex plane that encircles the origin a:: = in 
such a way that all singularities of F{£^, x) (as a function of x for fixed ^) are 
outside the contour. 

From Eq. (|54|) and the definition of F, we deduce that F(^, a;) satisfies the 
following partial differential equation: 

1 



2.9F ^dF 



X 



2^- 1 



F (57) 



This equation together with the requirement that x) ~ x + O(x^) as a; — > 
(as Gi = 1) determines uniquely F{£^,x). We find: 

(1 + x)^/l — x'^ 

G„ [C) is given by the contour integral in Eq. ([55| where the contour C encircles 
a; = in such a way that a: = 1 and x ~ ^\ are outside of the contour. 

Let us compute Gn{C) with ^ = Naz^ for large N , fixed z and n = cN with 
fixed < c < 1. We have 

G„(0 = :?^^^^(C,a;) - i ^ ; (59) 

where 

^ , , 2Q:z^a; , . , n 

$c(a;) = c Inx with c = — (60) 

is of order one for large N when n = cN with c of order one. For fixed c < 1 
and for large N , the contour integral can thus be computed using a saddle point 
method. The integral will be dominated by the neighbourhood of x* such that: 



dx 



„ . 2az n c 

= ie -t: = = — (61) 

{l + x*Y Nx* X* ^ ' 
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There exists a real solution for x* iff > — ie > . We want this condition 
to be satisfied for all n < N, therefore we must have > — , i.e. z > \ —. Our 



method can only describe the regime z > 



iey > 



which corresponds 



to the right tail of the distribution P(A,nax < y) (region where A,nax is above its 



mean value). Let us call the critical point j/cr = y —, i-e 



For y > ycr, there are two real solutions x* = —1 



Naz^ 



2n 
Naz^ 



The contour C must encircle but not 1 and —1, therefore we impose — 1 < 
a;* < 1. This implies 



Naz^ 



2n 



2n 

Naz^ 

2n 



Naz^ 
2n 



1 - 



2n 



2n 
Naz"^ 



(62) 



Thus 



az 



2c 



2c In 




-rrr-xHr-^} (63) 



The saddle point gives for large N: 



27r {l + x*)Vl-x* 



\ 



2tt 



N 



(64) 



where 



dx-^ 



In this subsection we have found, as given in Eq. (|64| . the expression of Gn(?) 



and thus the solution c„(?/, a) 



GniO (with ^ = Naz^ = ay^) of the 



recursion relation (ISU]) for large N and n — cN with fixed < c < 1. We have 
also shown that the validity of our approximation is the regime y > ycr with 



4.3 Computation of the distribution of Amax for large 
We want to compute for large N and for y > ycr (with the scaling y — 



large N) the cdf Pjv (A^^x < y) 



ZN{y,a) 
2'iv(oo,a) ' 



N for 

Using Eq. and the definition 
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of Gn in Eq. ([55)) . we get: 



InPw (Ainax < y) = lnZAr(y, a) - lnZAr(oo,Q!) 



2V< 



1 + 2^E 



GniO (65) 



n=l 

Therefore we need to compute In{0 = Y^n=i i^^lT^) ^"(C) for large N and 
^ = Naz^ (with fixed a and z). For that purpose, we do not actually need to 
use the approximate expression of G„ for n of order TV that we derived in the 
previous subsection. We can use the formal expression of G„ as a contour inte- 
gral G„(^) = §^ ^ ^J-^i F{£^^ x) and compute the sum over n before computing 
the integral by saddle point method. In particular we have: 



N-l 



dx F{i,x) 



dx F{^,x) 
7n 



c2iTrx{x-l) Jc 2iT: x'^ (1 - x) 2i7r - 



dx F{^, x) 



The function 



-t — TT — — , e="+i has indeed no singularity at the origin, 
its integral is thus zero. On the other hand, we have: 



E 



G„(e) 



dx F{i,x) 2Fi(l,l-iV,2- Ar,a;) 



2i7r X 



N 



N - 1 



where 2^1 is a hypergeometric function: 2Fi{a,b,c, z) = ^ 
(a)fe = a{a + l)...(a + fc — 1). For large N, we find: 



{a)k{b)k -^-i 



TV- 1 



fc>0 



k>Q 



N 



1 



+ 



1 



iV(l-x) (l-a;)2iV2 



Therefore we get for large N 



I f dx F{^,x) 



(66) 



Equivalently we can write: 



1 



dx 



7V/c 2tTT{l + x)i{l-x)^ 



(67) 



where $i(a;) = '^"^i ^ In a; (see (pU)) ). For large N, the saddle point method 
thus gives 



1 1 



N2tt (i + x*)i{l- x*y. 



\ 



2tt 



N 



(68) 
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where x* is given in Eq. (|^^ with n = N: 



Thus we find 



-1 + az^ 



1 - 1 1 



1 1 



az 



1 



TV 2 ^/2^^ AVaz^iaz"^ - 2)5 



1 - 1 1 



(69) 



(70) 



with 



$i(a;*) = az"^ 



l-\ 1 



az'- 



-21n (71) 



Therefore 



hlPAT (Ainax < 2/) = hiZ7v(2/, a) - lnZAr(cxD,a) 



{1 + 2^/^(0} 
1 - 



7V2V2^(a22 - 2)i 



As $i(x*) > for z > Zcr, i-e. olz^ > 2, the first term in the parenthesis can be 
neglected for large N: 



InPjv (A„,ax < y) 



-Nc 



1 



47r7V ^/2(az2_2)i 



(72) 



with y — zyN. Therefore we get the expression of the right tail of the cdf of 



-i e-^^'^+W for z>J- 

47rA^V2(az2 - 2)5 Va 



where the rate function ip+i^) = — — f^^^ ^ is given by: 



az^ 

V'+(z) = — 





+ ln| 




V az^ 





- 1 



(73) 



(74) 



We have thus found 

(Amax < zVn) ~ 1 



N 



AirNV^iaz^ - 2)' 



> 



(75) 
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Deriving with respect io t — z\/ N we get an equivalent of the probabihty density 
function of Amax for large N: 

V (A^a. = « e-^^^- (^) 

^ ' 2TTV2{at^ -2N) 

(76) 

We thus recover the dominant order for large N (given by V-'+(-2^)) that was 
derived in [5F by a Coloumb gas method, but in addition this method also 
provides us with the first correction term to the dominant order. 
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Let us now see how this precise right tail large deviation behavior in (|75p 
behaves when z from the right. Using the leading expansion of ip-i-{z) 

around z = \/2 in (HT]) and setting y = -y/^ + iV~i/6_|^^ i ^., z = y/^N = 

1 ^ ]^-2/3x Qj^g gjjjg f^Qj^^ ^ 



/2q 



'n (A,„ax <y)~l- — ^ e-l-^ . (77) 



_ 2 3/2 
3^ 



On the other hand, using the boundary condition q{z) « Ai(2) « ^^^^^y^ e 
as z — oo in the definition of the Tracy- Widom function one can easily de- 
rive the precise leading asymptotics of its right tail, F2(x) w 1 e~3^' 
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as x — >■ oo. Thus, our right large deviation function for small argument (when 
the fluctuation of A^ax to the right of its mean value \J2N ja is of 0(iV~^/^)) 
in ((77|) . matches smoothly with the precise right tail of the Tracy- Widom dis- 
tribution F2{x). 



5 Double scaling limit and Tracy- Widom distri- 
bution 

In this section, we provide an elementary derivation of the Tracy- Widom law 
for the GUE based on simple scaling analysis of the recursion relations derived 

in Section 2 in the vicinity of the critical point y — ya = \f^- This derivation, 
in our opinion, is mathematically simpler than the original derivation by Tracy 
and Widom [6l as it avoids the sophisticated asymptotic analysis of Fred- 
holm determinants. The derivation of the Painleve II equation from the scaling 
analysis of recursion relations that we follow here is similar in spirit (though 
rather different in details) to the analysis of the partition function in the two 
dimensional Yang-Mills theory on a sphere by Gross and Matytsin 391 . 
Let us recall that the Tracy- Widom distribution -^2(2^) is defined as 

F2{x) = exp | - y ds[s - x)g2(s)| (78) 
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where q{x) satisfies the Painleve II equation with the boundary condition 
q"{x) = 2q^{x)+xq{x) Painleve II 
'^^'^^ as x^oo (79) 

From Eq. it follows that '^^^'^^"'^ = -q^ix). 

We want to show that for large N the probabality of small typical fluctua- 
tions of Amax around its mean value are described by the Tracy- Widom 
distribution. For this, we need to first estimate how do these typical fiuctuations 
scale with N for large N. In the vicinity of the mean ^ let us then write 

(80) 

where is the scale of the typical fiuctuation and the random variable x has 
an N independent distribution for large N. Evidently the exponent 7 < 1/2 
(so that the fiuctuation is less than the mean) whose precise value is yet to 
be determined. Note also that since Amax always appears in the distribution 
IPAr(Amax < y) in the scaling combination yy/a (see Eq. pO)) ). we have chosen 
the prefactor of the fluctuation term as 1/ \/2a which then ensures that the 
random variable x describing the typical fluctuation is also independent of a. 

One way to estimate the exponent 7 is from the right large deviation tail 
computed in ([75t in the previous section. The right tail in (|75p describes the 
probability of large fluctuations of 0{\/N) to the right of the mean. Assuming 
that the right tail behaviour continues to hold even for fluctuations smaller than 

//V, we substitute z//V = + -7= iV^ a; in ([75]) . This gives 




N 



g^ + -LiV^.Ul- \ e-|-----^-^ (81) 



valid for x > 0, x large. Assuming that this continues to hold even for not 
so large x (so that it even captures the tail of the distribution of typical small 
fluctuations), we expect that in terms of this rescaled variable x, the tail of 
the distribution in (|81]) is independent of N for large N . Clearly, for this to 
happen the power of N must be zero both inside the exponential as well as in 
the prefactor in ((8T|) . indicating that \ + ^ = 0, thus 7 = — i. Hence, the 
correct scaling describing typical fluctuations, for large iV, is given by 

-^x (82) 

where x has an N independent distribution that we now have to compute and 
show that it is given by the Tracy- Widom function F2{x). 
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The meaning of double scaling limit is now clear. It simply says the follow- 
ing. Consider the cdf Pat (Amax < y) or rather its logarithm (for convenience) 
In Pat (Amax < J/)- In general, it is a function of two variables y and N. How- 
ever, in the vicinity of the mean y — > \J^^, if one takes the limit y — \J^^ — > 

and N ^ oo, but keeping the scaling combination x — \/2aN^^^ (y — ^) 
fixed, this function of two variables collapses into a function of the single scaled 
variable x 

In Pw (Amax < y) ^ ./ ^ N^'^ ^ (83) 

and our job is to to show that this scaling function f{x) = lni^2(a;) where -^2(2^) 
is the Tracy- Widom function defined in ([78]). In other words, we want to show 
that f"{x) — —q^{x) where q{x) satisfies the Painleve II equation ([79l) . 

Our starting point is the definition of the cdf in From Eq. it is 

easy to see that the partition function Z^r satisfies the recursion 

ZN^i{y,a) ZN+i{y,a) hN{y,a) 

^27 N = T 7 7 = RN[y, a) (84) 

Taking logarithm and using the definition in (|28p we get 

InP^v+l (Amax < y)+lnPAr_i (Amax < y)-21nPAr (Amax < 2/) = In f ^^^^ 

(85) 

In the double scaling limit, we will now substitute the anticipated scaling form 
in ((83| for the logarithm of the cdf on the left hand side of (|85|) . But we need to 
first evaluate lnPAr±i (Amax < y)- Replacing TV by TV ± 1 in (155)) and expanding 

for large TV, with x — \/2a N^/^ (^y — \J fixed, we get 



InPjvii (Amax < 2/) = f^^V2^{N±lf''^-^j'^^^^^ 

= f {xtN-^'^± — ± - + . 

\ QN 12 

= f{x) T N-^'^J'ix) + ^— /"(x) + 0(TV-i)(86) 



Substituting this result in (|85|) we get for the left hand side 

InPjv+l (Amax <y)+ lnPjV-1 (Amax < V) - 2 In P^ (Amax < 2/) 

kN-^'^ f"{x)+0{N~^) (87) 
From Eq. and dH?]), we get for large TV 

TV-^/3 fix) « In ( 3E^\ « In ( ^E^\ (gg) 
^ ^ ' \RNioo,a)J \ N/{2a) ) ^ ' 
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as i?Ar(oo,a) = N/{2a) (see Eq. dS])). This suggests that in this scahng hmit, 
Rn must scale as Rn {y, Q^) ~ (l + f"{x) + ...^ . More precisely, this 
leads us to the following large N expansion of RN(y,a) in the double scaling 
limit 

Rn {y, «) « ^ f 1 + N--^ n{x) + N-\2{x) + iV-^r3(a;) + ...) , (89) 
where 

f"{x)=ri{x) (90) 

and 7'2(a;), r^{x) etc. describing the higher order scaling corrections. Thus, if we 
can now determine the first subleading scaling function r-^ix) in the expansion 
of Rn iy,cx), then we can determine f{x) by integrating ri{x) twice. So, our 
next task is to determine ri{x) by analysing the recursion relations (|4T|) and 
(|42|) (setting n = N) in the double scaling limit. 

We now know, from (j89p . how i?jv(2/, Qf) behaves in the scaling limit with the 

scaling combination x = \f2a N^/^ (jj — fixed. In order to analyse the 

recursion relations (|4ip and we also need to know how SN{y, ct) behaves in 
this scaling limit. In order to match the leading N behavior of R]\[{y, a) with x 
fixed in (|42|) . it is not difficult to see that to leading order for large N, S^iy, a) 
must have the following scaling behaviour 

SN{y,a) « ( V2^N^" iy~\f^\\ +0{N-^'^) (91) 



2a ^ 

where si{x) is the leading order scaling function. Let us first evaluate the 
difference SN~i{y,o) — SN{y,a) that appears in (H^ . Replacing iV by iV — 1 
in (|9T|) , expanding for large N, we get 

/V-l/2 

SN-i{y,a) - SN{y,a) « —^s[{x) + 0{N-^/'^). (92) 

V 2a 



It rests to evaluate the partial derivative in (H^ . From the definition 

of the scaling variable x = \/2aN^^^ (^y — ^ , it follows, using chain rule, 

91ni?Ar(j/, a) d\nRN{y,a)dx 
dy dx dy 

= ^^i/e^lniMy^ 
dx 

= V2^N~^^^r[ix) + 0{N^^/^) (93) 

Finally substituting and in (with n = N) we get, 

V2^ (r[{x)N'^ + O(iV-t)) « = 2a(5jv-i-5jv) « V2^ (a^"^s'i(x) + ©(A^^^))] 
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Matching the leading order N~'^/'^ term gives a relation between si{x) and ri {x) : 
Si{x) = r[{x), i.e., 51(2:) = ri{x) + cq with cq a constant. From p4|) and the 
fact that when y 00, x 00, it follows that both the scaling functions ri{x) 
and si{x) must vanish as a; ^ 00. Thus the constant cq = and we have, for 
all X, 

ri{x) = si{x). (94) 

Having determined the relation ri{x) ~ si{x), we need one more relation 
between these two scaling functions in order to determine them individually. 
This will now be done by substituting the scaling solutions for i?jv(2/, ck) (given 
in ^89\i ) and SN{y, a) (given in (|9T|)') into the remaining recursion relation (|4T|) 



To analyse (HT|) (seeting n — N) , we need to evaluate the derivative ^'"''g^^^'"'' ■ 
From the definition of the scaling variable, x = \/2a N^^^ (^y — ^ , it fol- 
lows, that = ^ + ■ We then use the chain rule and to express 

91n_Rjv(w, a) r\ (x) — 1 N^^ , , , A^^i ^ 1 , „ / n / . n „ / . m 
= + ^20^ +^ b^'^ a: - 2ri X r X + 2r^ x +... 

Again replacing iV by ± 1 in (|89l) and expanding for large A^, keeping x = 
V2aN'^/^ [y - fixed, we get 

r' ix) — 1 iV^^ A^~i 

Rn-i-Rn+i - + r'2{x) + —- [-2ri(x) - xr[{x) + 6r!,(a;) + r'/'(x)]+.. 

a a ba 



and similarly from (|9ip 

-S'at-i - "S*^ ~ ^^^^ — si(x)si(x) + ... 
a 

Substituting these results in (|4T|) and matching the leading order term {0{N^^/'^)), 
we get the desired second relation between ri{x) and si{x) 

xr[{x) - 2n{x)r[{x) = -^n{x) ~ ^r[{x) + ^r'l' {x) + 2si{x)s[{x). 

Eliminating si(x) by using si(x) — ri{x) we finally get a single closed equation 
for ri {x) 

2xr[ (x) + ri (x) = ^r'/' (x) + 6ri (x)ri (x) . (95) 

Let us write 

ri(x) = -M^(x) (96) 
Eq. (j95p then becomes an equation for u{x): 

u{u"' - 6u^u' - xu' -u) ^ ~3u'{u" -xu- 2v?) (97) 

Let W{x) = m"(x) - xu(x) - 2u^{x). Then dST]) becomes 

^ dWXx2^_g^ (98) 
ax 
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which can simply be integrated to give 

A 

u{x) 



W{x) = -^j-^. (99) 



where A is an arbitrary constant. Hence we have 

u"{x) - xu{x) - 2u^{x) = — (100) 

u{xy 

From the boundary condition ri(x) — as a; — ^ oo (which follows from (01])), 
it follows using ri(x) = v?{x) that u{x) — as a; — )■ oo. Taking x — oo in 
(jlOOp then fixes the value of the constant ^ = 0. Finally, from (|90)) . we have 
f"{x) = ri{x) = — tt^(a;) where u{x) satisfies the Painleve II equation 

u"{x) ^ xu{x) + 2u^{x) (101) 

To fix the boundary condition for u{x), we again invoke the matching with 
the right large deviation tail in (|8T|) . Taking logarithm of (|8T|) with 7 = —1/6 
and using lnPjv(Ainax ^ y,a) — f{x) we find that 

/(x)« ^-e"^^^ as a;^cx) (102) 

167ra;2 

Hence u^(x) = —f"{x) e~3^^ i^TT? ^^"^ consequently as x — ?► 00 

«(.) = V^m«e-i-^^^^. (103) 

Finally integrating f"{x) twice and using the appropriate boundary condi- 
tion as X ^ 00, we get 

POO 

f{x)^- ds{s~x)u^{s) (104) 

where u{x) satisfies the Painleve II equation (|10ip with the boundary condi- 
tion (|103p . Comparing to ((75)) . we have thus shown that the scaling func- 
tion f{x) — lni^2(a;) where ^2(2;) is the Tracy-Widom function (/3 — 2). This 
then constitutes our derivation for the Tracy-Widom distribution for the GUE 
{13 = 2). 



Using recursion relations that we have derived for orthogonal polynomials 
on a semi-infinite interval, we have shown that the large N asymptotics of the 
distribution of the maximal eigenvalue of a Gaussian random matrix (from the 
GUE) is described in the double scaling regime by the Painleve II equation. 
Similar recursion relations for other orthogonal polynomials leading to different 
kinds of Painleve equations have also been established in a number of papers, 
see [15] and references therein, in particular [SO] HE] HZ] . 
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6 Conclusion 



In this paper, we have provided a rather simple and pedestrian derivation of the 
Tracy- Widom law for the distribution of the largest eigenvalue of a Gaussian 
unitary random matrix. This was done by suitably adapting a method of orthog- 
onal polynomials developed by Gross and Matytsin [39 in the context of two 
dimensional Yang-Mills theory. Our derivation requires just elementary asymp- 
totic scaling analysis of a pair of coupled nonlinear recursion relations. Strictly 
in the N ^ oo limit, there is a 3-rd order phase transition in the form of the 
probability distribution of Amax as Amax crosses its mean value from left to right. 
For finite but large the two regions are connected by a smooth crossover 
function and the shape of this crossover function is precisely the Tracy- Widom 
distribution that describes the 'typical' small fiuctuations of Amax around its 
mean. The 'atypical' large fluctuations to the left and right of the mean are 
described by large deviation tails that correspond to the two 'phases' across 
this phase transition. In qualitative analogy to the two-dimensional Yang-Mills 
theory, the left (left large deviation) and the right (right large deviation) phases 
correspond respectively to the 'strong' and 'weak' coupling phases of the two- 
dimensional QCD. Apart from the simple derivation of the Tracy- Widom GUE 
law, we were also able to compute the precise right large deviation tail of the 
maximal eigenvalue distribution that is not described by the Tracy- Widom dis- 
tribution. In the language of QCD, this right tail corrections are similar to the 
non-perturbative (in expansion) corrections to the QCD partition function 
in 2-d I3S]. 

One drawback of our method is that it works only for the GUE case (with 
Dyson index P — 2). It would be challenging to see if this method can be ex- 
tended/generalized to derive the Tracy- Widom law for the other two Gaussian 
ensembles, namely the GOE (/3 = 1) and the GSE {(3 — 4:). In addition, this 
method for /3 = 2 should be useful to compute the largest eigenvalue distribu- 
tion for other non-Gaussian matrix ensembles, such as the Laguerre (Wishart 
matrices) or the Jacobi ensembles. 
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